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The application of an efficient flow control system on helicopter rotor blades may lead to improved aerodynamic 
performance. Recently, our invention of Rod Vortex Generators (RVGs) has been analyzed for helicopter rotor 
blades in hover with success. As a step forward, the study has been extended to forward flight conditions. For this 
reason, a validation of the numerical modelling for a reference helicopter rotor (without flow control) is needed. 
The article presents a study of the flow-field of the AH-1G helicopter rotor in low-, medium- and high-speed for- 
ward flight. The CFD code FLOWer from DLR has proven to be a suitable tool for the aerodynamic analysis of 
the two-bladed rotor without any artificial wake modelling. It solves the URANS equations with LEA (Linear 
Explicit Algebraic stress) k-@ model using the chimera overlapping grids technique. Validation of the numerical 
model uses comparison with the detailed flight test data gathered by Cross J. L. and Watts M. E. during the Tip 
Aerodynamics and Acoustics Test (TAAT) conducted at NASA in 1981. Satisfactory agreements for all speed re- 
gimes and a presence of significant flow separation in high-speed forward flight suggest a possible benefit from 
the future implementation of RVGs. The numerical results based on the URANS approach are presented not only 


for a popular, low-speed case commonly used in rotorcraft community for CFD codes validation but preferably 


for medium- and high-speed test conditions that have not been published to date. 
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Introduction 


CFD simulation of a rotorcraft is very challenging 
mainly due to the unsteady phenomena occurring during 
a revolution (e.g. compressibility effects, shock waves, 
dynamic stall and blade-vortex interactions). Thus, an 
accurate simulation of the flow-field passing a helicopter 
rotor blade is essential for proper prediction of exerted 
airloads. The output of the Tip Aerodynamics and Acous- 
tics Test (TAAT) conducted by Cross and Watts [1] in 
1981 for the AH-1G helicopter is very popular within 
rotorcraft community for CFD codes validation. It com- 
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bines an extended set of flight test data, such as overall 
thrust and torque coefficients, instantaneous values of 
normal force and pitching moment coefficients at several 
cross-sections as well as pressure distributions for se- 
lected span-wise locations and azimuthal positions of the 
blade. A literature survey of the published numerical re- 
sults [2-8] for the TAAT AH-1G campaign reveals that 
the low-speed flight (V.. = 150 km/h) has been simulated 
in the past by several researchers with success, in con- 
trast to the higher speed flight test points where CFD 
calculations are very limited (i.e. medium- and high- 
speed flights with forward velocities of V.. = 210 km/h 
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Nomenclature 


AR rotor aspect ratio [-] 


CFL Courant Friedrichs Lewy number [-] 
Cn pitching moment coefficient [-] 

Ca normal force coefficient [-] 

Cp pressure coefficient [-] 

Cr, Cg thrust and torque coefficients [-] 

c blade chord [m] 

Mr, Ma tip and forward flight Mach number [-] 
Pr Prandtl number [-] 

P, static pressure [Pa] 

Po atmospheric pressure [Pa] 

R rotor radius [m] 

Rgas gas constant [J/kg/K] 

r radial position [m] 

Rer tip Reynolds number [-] 


and 290 km/h respectively). Various numerical ap- 
proaches/methods have been investigated over the last 
years concerning the low-speed case. Hernandez and John- 
son [2] used the FPR-CAMRAD/JA code supplemented by 
a rotor wake model, while Ramachandran et al. [3] de- 
veloped a method based upon the vorticity-embedding 
technique. Ahmad and Duque [4] applied the chimera 
method to account for the blade motion in their research. 
More recent CFD calculations can be found in [5-8]. 
Hernandez et. al. [2] and Tan et. al. [8] presented unique 
numerical results for the high-speed forward flight case 
based on an inviscid approach supplied with external 
wake model. In contrast, our method does not require any 
prior knowledge of the flow-field. It is sufficient to 
supply the flight conditions (i.e. angular and forward 
flight velocities and blade kinematics) as input. For this 
reason, the paper presents numerical results based on the 
URANS approach for the medium- and high-speed cases, 
which are unique in terms of the literature survey. As a 
first approximation the helicopter rotor blades are iso- 
lated and considered rigid. Still, the remaining task is 
computationally very demanding. 

In the past, the authors of the paper have been dealing 
with numerical simulations of helicopter rotor blades in 
hover and forward flight conditions for various rotors 
(for example Caradonna-Tung [9-10] or PZL W-3A Fal- 
con (Sokol) [11]). The initial results of the simulations of 
the AH-1G rotor have been previously published in [12]. 


Geometrical Model and Chimera set-up for an Iso- 
lated Rotor of the AH-1G Helicopter 


The AH-1G is a two person, single-engine helicopter 
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Vo forward flight velocity [km/h] 
x chordwise location [m] 
X,Y,Z Cartesian coordinates [m] 
Greek symbols 
B flapping angle [°] 
Bo Bs backward and sideway disk tilt [°] 
Bo coning [°] 
y specific heat ratio [-] 
83 rotor pitch-flap coupling [°] 
8 pitching angle [°] 
ba 9; lateral and longitudinal cyclic coeff. [°] 
Oo, 8.75 collective at blade root and r/R = .75 [°] 
u advance ratio [-] 
o rotor solidity [-] 
y blade azimuthal position [°] 


that first flew in 1965 and has been out of service for 20 
years. The AH-1G 2-bladed, teetering rotor with a 540 
symmetrical airfoil section, which is a highly modified 
derivative of the NACA 0012 (thickness ratio = .0933), 
has a rectangular planform with a chord c = 0.686 m. The 
blades have a linear twist of -10° from the shaft to the tip 
(effective twist from blade root to tip is —8.48°) with a 
radius R = 6.71 m (aspect ratio AR = 9.8). During the 
TAAT test, the effective chord of the blade was increased 
to accommodate the instrumentation sleeve, reducing the 
effective aspect ratio to 9.2. 

In forward flight, the blades not only rotate around the 
azimuth y (as for hover conditions), but also pitch (0(y)) 
and flap (A(y)). The azimuthal angle of 0° refers to the 
position of the first blade pointing in the direction oppo- 
site to the flight velocity (therefore y = 0° — 180° denotes 
the advancing side while y = 180° — 360° marks the re- 
treating side). The recorded in-flight blade motion is de- 
scribed in a form of Fourier series using first harmonics 
which are tabulated in [1] and applied in the simulations: 

0 = bo + O. > cos (y) + 0, > sin (y) (1) 
B = Bo + Be > cos (y) + Bs > sin (y) (2) 

The computational domain is formed by three chimera 
overlapping component grids. This technique allows for 
the easy generation of grids for complex configurations 
by decomposing them into simple, independent parts. 
Using this approach for our study, it is relatively easy to 
implement the blade motion without any grid deforma- 
tion and to add chimera component grids designed for 
flow control devices (i.e. RVGs already analyzed for ho- 
vering rotors in [10]). The component overlapped meshes 
for the AH-1G helicopter rotor has been generated semi- 
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automatically using a python based internal programming 
language embedded within the Interactive Grid Generator 
(IGG) — a commercial software from Numeca Interna- 
tional. A component grid is created for a single rotor 
blade placed within the background mesh (the remaining 
blade is generated automatically). The Cartesian back- 
ground component grid (Fig.1) is designed with dimen- 
sions of 16.4 R x 18.2 R x 18.2 R, consequently the outer 
part of the domain is located at least 8.0 R away from the 
rotor disk in every direction. The vicinity of the rotor and 
its wake is solved with a cuboid structure with uniform 
volumes (0.1 c x 0.1 c x 0.1 c) increasing the accuracy of 
interpolation between blade component and background 
grid, as well as support acoustic wave near-field propa- 
gation up to a frequency of approximately 1 kHz. Alto- 
gether, 32 computational blocks contain 9.4 x 10° control 
volumes. 

The blade component grid (Fig. 2) is of C-type in 
streamwise and H-type in spanwise direction. It spans 
from the surface 1.2 c in all directions (normal and radi- 
al). Designated 40 computational blocks contain 3.9 x 
10° of volumes per blade. The non-dimensional distance 
of the first layer of cells from the solid surface of the 
blade is of the order of 1 which is sufficient for resolving 
the turbulent boundary layer using low-Re turbulence 
model of LEA k-œ down to the wall. Figure 3 presents a 
combined chimera overlapping grid topology for the 
AH-1G helicopter rotor. Overall, the whole chimera set- 
up is formed by 3 component grids (background and two 
blades) with 112 blocks and 17.2 x 10° volumes. 


Fig. 1 Background chimera component grid 


Three types of boundary conditions are applied in the 
numerical simulation: no-slip condition with zero heat- 
flux (adiabatic) at the surface of the rotor blades, far-field 
condition at the exterior of the background grid and a 
special chimera condition at the outer edge of the blade 
component grid which is necessary for the interpolation 
of flow variables between meshes. 


i Blade chimera component grid — >H — 


A : TK 


Fig.3 AH-1G rotor chimera overlapping grid topology. 


Physical and Numerical Modelling 


The present investigation was carried out with the 
FLOWer solver from DLR [13]. It is a modern, parallel, 
block-structured, cell-centered, aviation oriented CFD 
code solving the Favre-averaged Navier-Stokes equations 
with various turbulence models. The two-equation, low- 
Reynolds number, Linear Explicit Algebraic stress model 
(LEA k-w) was used for the numerical investigation due 
to improved transonic flow predictions as well as lack of 
the wall distance calculation in the closure (computation- 
ally demanding for the proposed approach). For the equ- 
ation of state, the perfect gas law is used with the specific 
heat ratio y = 1.4, gas constant Roas = 287.1 and Prandtl 
number Pr = 0.72. The laminar viscosity is calculated 
according to the Sutherland’s law. 

The numerical algorithm is based on a semi-discrete 
approach with a finite-volume, central scheme of 2”? or- 
der of accuracy for the spatial discretization. An explicit 
Runge-Kutta method (CFL = 2.5) of time integration was 
applied in the internal iterations of the time-accurate im- 
plicit dual-time-stepping scheme of 2" order of accuracy. 
The scalar artificial dissipation model of Jameson is im- 
plemented to damp numerical oscillations. A time step 
equal to 0.25° of azimuth (1440 time steps per revolution) 
and a drop of density residual by 3.0 orders of magnitude 
proved to be sufficient to obtain an accurate, unsteady 
flow-field around the rotor. 
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TAAT Forward Flight Cases Simulated 


After the TAAT flight test campaign, one hover IGE 
(test point no. 2370) and six forward flight cases (test 
points no. 2152-2157) with increasing forward speeds 
(from 150 km/h to 290 km/h) were documented in [1]. 
This article presents exemplary numerical results ob- 
tained for a subset of 3 test points (see Table 1). Test 
point no. 2157 (low-speed case, 150 km/h) was computed 
and validated against flight test data and compared with 
other CFD results found in the literature [2-8]. Two high- 
er speed test cases investigated do not have similar cita- 
tions available in the literature. Numerical results for test 
point no. 2152 (high-speed, 290 km/h) are satisfactory in 
terms of flight test data comparison revealing significant 
flow separation (shock-wave induced at the advancing 
side and due to dynamic stall at the retreating side). Mi- 
tigating the effect of this type of flow separation is a goal 
of our invention of Rod Vortex Generators (RVGs) [14]. 

Blade pitch is one of the key aircraft measurements 
needed to produce accurate CFD predictions. Given that 
the blade pitch angle varies due to the linear twist of the 
rotor blade of —10°, the location of the angle measure- 
ment is critical. During the TAAT, a sensor mounted be- 
tween the fixed hub and the feathering bearing housing 
measured blade pitch. CFD predictions generally use 
blade pitch values at 0.75 R, while the measured values 
from TAAT are relative to the blade root at 0.15 R (see 
Table 1). 


Table 1 TAAT flight test conditions [1] 


low-speed medium-speed high-speed 
point 2157 point 2155 point 2152 
Va 150 km/h 210 km/h 290 km/h 
M» 0.12 0.17 0.24 
Mr 0.65 0.65 0.64 
Rer 9.72*10° 9.78 10° 10.22x10° 
u 0.19 0.26 0.38 
80 11.7° 13.4° 18.0° 
0.75 5.7° 7.4° 12.0° 
0. 1.7° 2.12 3.6° 
ðs -5.5° -7.9° -11.8° 
Be 2.13° 2.38° 1.13° 
Bs -0.15° 0.19° 1.11° 
Cr 0.00464 0.00464 0.00474 


Conversion of blade pitch from the root measurement 
to 0.75 R is obtained with the following equation: 0.75 = 
o — (.75 —.15):(blade twist per meter radius)-R, where 
blade twist per meter is 1.493. Therefore, the values of 0o 
transferred to 0.75 R are equal to 5.7°, 7.4° and 12.0° for 
the low-, medium- and high-speed markers, respectively. 

In the course of conducting this study, great effort was 
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expended to ensure that the results are as accurate as 
feasible. The AH-1G rotor design includes a little known 
feature, called pitch-cone coupling, which reduces blade 
pitch as the rotor coning increases. This feature is not 
often discussed or included in analysis of this rotor. It 
was studied here to ascertain what impact this feature 
could have on blade pitch over these flight conditions. 

Pitch-cone coupling plays a role in keeping the rotor 
from stalling out as the helicopter pulls high g-loads. 
This teetering hub has a built in 2.75° precone to minim- 
ize flexing of the hub during normal flight. The hub de- 
sign incorporates a thin flat yoke, with low flapwise 
stiffness, just outboard of the flapping hinge located at 
the shaft. When the rotor experiences high loading (e.g. 
pulling g's in a maneuver), this region flexes, much as a 
rigid rotor flaps. Because this flex area is not in align- 
ment with the pitch rods, a de facto 53 angle (rotor pitch- 
flap coupling) is created. It should be noted here that the 
AH-1G teetering hub has no pitch-flap coupling, which is 
the usual result of hubs that include 63. For inelastic 
blade analysis, such as undertaken in our report, the pub- 
lished pitch-cone coupling ratio term for the AH-1G is 
equal to -0.682 [15]. 

As the test conditions studied here are lightly loaded 
(Cy/o = 0.068), the rotor forces could reduce rotor coning 
from the built-in precone. Basic calculations show that an 
articulated rotor, at this Cr /o, would cone in the range of 
2.1°. Due to the yoke stiffness in this teetering hub, def- 
lections from the precone angle will be smaller than the 
predicted coning reduction of 0.65°. Both elastic and 
inelastic rotor analysis of an AH-1G in level flight at 214 
km/h, were conducted in 1980 [16]. The inelastic model- 
ing predicted a blade pitch change of 0.03° while the 
elastic modeling predicted 0.55° due to pitch-cone cou- 
pling. In either case, higher speeds will result in coning 
closer to the precone 2.75°, resulting in reduced blade 
pitch deflections. Since the current study includes inelas- 
tic blade modeling, the limited elastic modeling coupling 
result was not used. Given that the inelastic coupling 
term produced such a minimal effect, pitch-cone coupl- 
ing was not included in the study. 

Due to relatively high uncertainty in collective reading 
based on the cabin lever position of the order of + 0.75° 
an additional verification has been conducted to assure 
zero normal force (zero pitch) at certain experimental 
azimuthal and spanwise position along the blade. Ac- 
cording to TAAT report (Table 1) the collective 89 correc- 
tion was set to be 0.3° for all cases under consideration 
being within bounds of the reading accuracy [1]. As a 
result, for the numerical simulations 8 7; was set to 6.0°, 
7.7° and 12.3° for the cases studied. 


Numerical Validation against the AH-1G Flight Test 
Data 


The thrust coefficient Cy is satisfactorily reproduced 
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by the CFD techniques employed here, when the flow 
conditions and blade motion reported in TAAT are ap- 
plied in the simulations. The thrust coefficient is only 
overpredicted by 2% for the low-speed case (Cr = 
0.00473). The thrust coefficient over-prediction increases 
to 9% (Cr = 0.00507) for the medium-speed test point. 
Finally, for the high-speed case the numerical simulation 
of the complex flow around the rotor increases the devia- 
tion in thrust up to 20% (Cr = 0.00569) [14]. Similar 
behavior was identified with analogous modelling and 
comparable flow conditions for the 4-bladed rotor of the 
PZL W-3A “Sokol” (Falcon) helicopter [11]. Initial re- 
view suggests reasons for the growing over-prediction 
trend are two folds: lack of the blade flexibility and rotor 
trimming in the modelling. 

As an example, Figure 4 (shown here using the Euro- 
pean customary clockwise rotation) presents the rotor 
wake (Q-criterion colored by vorticity magnitude) for the 
high-speed forward flight case. Although the wake is 
dissipated far from the rotor, the refinement levels are 
sufficient for properly reproducing the blade loading. The 
tip vortices and trailing edge vortex sheets are shed from 
the blades interacting in a complex manner. 


Fig.4 AH-1G helicopter rotor wake (high-speed flight) 


Figures 5—7 present exemplary numerical validation 
for all considered forward flight cases in terms of C, vs y 
and Cm vs y based on integration of pressure distribu- 
tions only (no viscous forces) at three cross-sections of 
the blade (r/R = 0.75, 0.96 and 0.99) which are repre- 
sentative of the flow-field inboards and outboards close 
to the tip. Figure 5 refers to the cross-section located at 
r/R = 0.75, the agreement with flight test data is satisfac- 
tory for the whole range of forward velocities analyzed. 
The maximum values of C, appearing at the retreating 
side are well predicted by CFD. Some detailed features 
of the flow presented in flight test data are well repro- 
duced by computations (for example local maxima of C, 
close to y = 90° for the low-speed case). 


Figure 6 presents the variation in time of normal force 
(C,) at span-wise station r/R = 0.99 close to the tip. As 
with the previous cross-section, the overall agreement 
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Fig.5 Instantaneous normal force coefficient C, at r/R = 0.75 
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Fig. 7 Instantaneous pitching moment coeff. Cm atr/R = 0.96 
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with flight test data is acceptable taking into account the 
number of simplifications in the modelling. The maxi- 
mum values of C, at the retreating side are well repro- 
duced by CFD. Close to zero normal force C, (the blade 
is pitching down close to the tip) caused by lack of inci- 
dence is present in the results of the numerical simula- 
tions as well. 

Usually, prediction of distribution of the momentum 
coefficient C,, is difficult for CFD even in cases of 2D 
airfoils. For this reason, and despite the number of sim- 
plifications applied to the numerical simulation presented 
in the paper, the agreement with flight test data is ac- 
ceptable for the entire speed sweep investigated at exem- 
plary span-wise station r/R = 0.96 (Fig. 7). The evolution 
of Cm is well reflected in the numerical simulations (the 
exception is the appearance of negative Cm close to y = 
0° for the high-speed case, not present in flight test data). 

It is worthy of mention that a correct prediction of 
integral parameters (like C, or C,,) do not ensure that all 
features of the flow-field have been well reproduced by 
the simulation (cancelation of errors), therefore during 
TAAT investigations, values of blade pressure coefficient 
C, were presented and analyzed. Figure 8 demonstrates 
example C, distributions (high-speed test point) for rep- 
resentative cross-sections and azimuthal positions. The 
shock wave existing at the suction side of the advancing 
blade (y = 90°) at r/R = 0.86 is captured by the numerical 
simulation in terms of pressure drop and shock location. 
For the radial station of r/R = 0.96, but at the retreating 
side (y = 270°), the comparison of C, distributions is also 
satisfactory. Lastly, two shock waves pattern present at 
the advancing side (y = 90°) for the tip section (r/R = 
0.99) is well reproduced by CFD. Both, flight test data 
and numerical results indicate non-lifting conditions due 
to the blade pitch down with a nearly zero incidence. 

The acoustic part of TAAT campaign included read- 
ings of instantaneous value of static pressure at several 
locations in accordance with the position of the pressure 
taps designed for capturing C, distributions. Figure 9 
presents a variation of pressure (P,/P,) at the blade suc- 
tion side at radial section of r/R = 0.86 and chordwise 
location of x/c = 0.25. A sudden drop of pressure in high- 
speed flight at the advancing side is caused by the dy- 
namic appearance of strong shock-wave also visible in 
previously described flight test data and CFD results of 
Fig. 8. When the shock-wave is alleviated, the pressure is 
recovered reaching a maximum at the retreating side. 

The disk loads in terms of normal (C,), torque (C,) 
and pitching moment (C,,) Mach scaled coefficients are 
presented in Figure 10. When increasing the forward 
flight velocity, more normal force is generated by fore 
and aft parts of the rotor disk decreasing contributions of 
the advancing and retreating sides. Advancing side nega- 
tive C, region due to nose-down pitch close to the tip and 
retreating side reverse flow area due to negative inflow 
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velocity grow with forward flight velocity. Higher trans- 
lational velocity of the rotor results in more pronounced 
areas of positive torque coefficient (C,) on the retreating 
side and negative C, close to the tip on the advancing 
side. In case of pitching moment coefficient (Cm), the 
majority of the positive influence is located on the aft 
part of the rotor disk shifting to the advancing side. For 
the high-speed forward flight, the numerical results pre- 
dict a negative value of C,, close to the tip at y ~ 0° that 
is not present in flight test data. 


r/R = 0.86, y = 90° 
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Fig. 8 Pressure coefficient C, distributions at selected radial 
sections and azimuthal positions for the high-speed 
forward flight case 
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Conclusions 


The paper presents the results of the numerical simu- 
lation of the isolated main rotor of the AH-1G helicopter 
in forward flight conditions reported in TAAT. The 
FLOWer code from DLR has been validated proving its 
capability for reproducing the flow-field of the rotor 
blades. The URANS results for the medium- and high- 
speed cases are unique in terms of the literature survey 
showing an acceptable agreement with the test data. The 
predicted existence of flow separation for the high-speed 
case reveals the possibility of future implementation of 
flow control systems (e.g. vortex generators [12], wall 
perforations [17] or synthetic jets [18]) aiming at im- 
provement of the aerodynamic performance. 
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Fig. 10 Mach-scaled normal force (C,), torque (C,) and pitching moment (Cm) coefficients in low-, medium- and high-speed for- 
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